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System-size dependence of the ground-state energy is investigated for A'^-site one- 
dimensional (ID) quantum systems with open boundary condition, where the interaction 
strength decreases towards the both ends of the system. For the spinless Fermions on the 
ID lattice we have considered, it is shown that the finite-size correction to the energy per 
site, which is defined as E^/N — linif^^^ E^/N, is of the order of 1/A^^ when the reduction 
factor of the interaction is expressed by a sinusoidal function. We discuss the origin of this 
fast convergence from the view point of the spherical geometry. 

§1. Introduction 

A purpose of numerical studies in condensed matter physics is to obtain bulk 
properties of systems in the thermodynamic limit. In principle numerical methods 
are applicable to systems with finite degrees of freedom, and therefore occasionally 
it is impossible to treat infinite system directly. A way of estimating the thermody- 
namic limit is to study finite-size systems, and subtract the finite-size corrections by 
means of extrapolation with respect to the system size.'^''^ 

As an example of extensive functions, which is essential for bulk properties, 
we consider the ground state energy of iV-site one-dimensional (ID) quantum 
systems. In this article we focus on the convergence of energy per site E^ /N with 
respect to the system size A'^. In order to clarify the discussion, we specify the form 
of lattice Hamiltonian 

H = Y,h,^,_,^ + Y,g,, (1-1) 

which contains on-site terms and nearest neighbor interactions h^i_^_^. We assume 

that the operator form of /i^ ^_|_^ and 5^ are independent of the site index £, which 

means that H is translationally invariant in the infinite N limit. It is possible to 
include g^ into h^ by the redefinition 

and therefore we group g^ with hg^_^^ as shown in Eq. (1-2) if it is convenient. A 
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typical example of such H is the spin Hamiltonian of the Heisenberg chain 

H = jJ2^i-^e+i-Bj2^f^ (1-3) 

e e 

where represents the spin operator at £-th site, and Sf its Z-component. The 
parameters J and B are, respectively, the neighboring interaction strength and the 
external magnetic field. In this case h^^^^ and are, respectively, JS^ • Sg_^i and 

—BSf. If the chain is infinitely long, H in Eq. (1-3) is translational invariant, and 
the ground state [iP'q) is uniform when there is no symmetry breaking. For example, 
the bond-energy J ( • ^^j^]) = J \ • S^^^^ I^'q) of the integer-spin Heisenberg 
chain is independent on I. 

This homogeneous property of the system is violated if only a part of the inter- 
actions h^2i ^2 3' • • •' ^"^^ N present, and the rest does not exist. In other 
words, if we consider an N-site open boundary system defined by the Hamiltonian 

7V-1 N 
Hopen = Yl hi+l + ^91, (1-4) 

the ground state I^^q) is normally non- uniform. As a result the expectation values 
(/i^^+i) and (g^) are position dependent, especially near the boundary of the system. 
The ground state energy of this A^-site system is normally not proportional to 
the system size N, which shows the presence of boundary energy correction. Such a 
finite-size effect is non-trivial when the system is gap less, as observed in the S = 1/2 
Heisenberg spin chain:-i^' 

In case that we are interested in the bulk property of the system, it is better to 
reduce the boundary effect as rapidly as possible. For this purpose Vekic and White 
introduced a sort of smoothing factor to the Hamiltonian 

N-l . g _^ g \ 

-^Smooth = X/ i^e,e+i^ 2 ) ' ^^'"^^ 

where A^ is almost unity deep inside the system and decays to zero near the both 
boundaries of the system.^J The factor A^ is adjusted so that the boundary effect 
disappears rapidly with respect to the distance from the boundary. A simplest 
parametrization is to reduce only A^ and ^^v-i unity, leaving other factors 

equal to unity. This simple choice of A^ is often used for calculations of the Haldane 
gap^ 

As an alternative approach, Ueda and Nishino recently introduced the hyperbolic 
deformation, which is characterized by the non-uniform Hamiltonian 

where A is a small positive constant of the order of 0.01 ~ O.l.l^'l^ As long 
as the form of the Hamiltonian is concerned, -f^Hyp ^^'^ regarded as a spe- 
cial case of -f^smooth -^l- (1'^) with A^, = cosh (A ^^~^~^ ) • But in the scheme 
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of hyperbolic deformation, the factor cosh (A ^— ^) is an increasing function of 
\(2i — N — l)/2\, and therefore the boundary effect is in principle enhanced. This 
enhancement works uniformly for most of the lattice sites, and the expectation value 
{h£i_^_i) = {^q \ I^^o) ground state |!?'o) becomes nearly independent on 

i for most of the bonds. After obtaining the expectation value {h^^_^_^) at the center 
of the system for several values of the deformation parameter A, one can perform an 
extrapolation towards A = to get the energy per site of the undeformed system. 
Such an extrapolation is possible since the hyperbolic deformation has an effect of 
decreasing the correlation length of the system. 

The hyperbolically deformed system is closely related to classical lattice models 
on the hyperbolic plane with a constant and negative curvature.'9''9''3''9'l3™'1D'II3'II3''2J'15J''lS) 
In this article we imagine the case of a positive constant curvature, where the clas- 
sical lattice models are on a sphere. The corresponding quantum Hamiltonian can 
be written as 

^Sph. = sm — ( /i^^^+i + 1 , (1-7) 



where A^^ = sm{i7r/N) decreases to zero toward the system boundary. We call such 
a modification of the bond strength as the spherical deformation, and consider N as 
the system size. We analyze the ground state \\Pq) and the ground-state energy 
of this deformed Hamiltonian for the case of spinless free Fermions on the lattice. 
We find that the difference 

lim , (1-8 



which is the finite-size correction included in the energy per site E^/N, is of the 
order of 1/N'^. Note that this dependence is the same as observed for the 

system with periodic boundary conditions, described by the Hamiltonian 

^Penod. = E [k^^l + ^^^) + (Vl + • (1-9) 



In a certain sense, the spherically deformed system does not contain system bound- 
ary. 

Structure of this article is as follows. In the next section we introduce a spinless 
free Fermion model on ID lattice. For tutorial purpose, the finite-size effect is 
reviewed for systems with open and periodic boundary conditions. In Sec. 3 we show 
our numerical results obtained from the diagonalization of the spherically deformed 
Hamiltonian -ffgpj^ in Eq. (1-7). In Sec. 4 we consider geometrical meaning of the 
spherical deformation by way of the Trotter decomposition applied to the deformed 
Hamiltonian. We also consider a continuous limit, where the lattice spacing becomes 
zero. We summarize the obtained results in the last section. 
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§2. Energy corrections in the free fermion system 

As an example of ID quantum systems, we consider the spinless free Fermions 
on the ID lattice. The Hamiltonian is defined as 

e e 

where t and /j, are, respectively, the hopping parameter and the chemical potential. 
For simplicity we set /x = and treat the half-filled state in this Section when is not 
explicitly shown. As a preparation for the spherical deformation, let us observe the 
ground state properties of the above Hamiltonian, when open or periodic boundary 
conditions are imposed for finite-size systems at half filling. 

First we consider the A^-site system with open boundary conditions, where the 
Hamiltonian is written as 

Af-l 

^o = -^ E {4^i+i + 4+i^e) ■ (2-2) 
e=i 

Since there is no interaction, the one-particle eigenstate |^^) represented by the 
wave function 

mel^J = = sin (2-3) 

is essential for the ground-state analysis, where m is the integer within the range 
1 < m < N. The corresponding one particle energy is 

rnvr 

= -2t cos , (2-4) 

and the ground-state energy at half filling is obtained by summing up all the negative 
eigenvalues. Assuming that N is even, the ground-state energy is obtained as 



N/2 

m=l 



tt/2 

sm ■ 



N + 1 



(2-5) 



after a short calculation. Expanding the r.h.s. with respect to A^, one finds the 
asymptotic form 

2 t /2 \ 
t + — (2-6) 



N TT N \Tr ^ 
Compared with the energy per site in the thermodynamic limit 

lim — ^ = - / -2t cos kdk = 1, (2-7) 

N^oo N TT Jo TT ^ ^ 

it is shown that the finite-size correction to the energy per site (or even to the energy 
per bond) is of the order of 1/N. 



Spherical Deformations for ID Quantum Systems 



5 



The A^-dependence of the energy correction changes if we impose the periodic 
boundary conditions, where the Hamiltonian is given by 

N-l 

Hp = -tY^ (cjc^+i + - * (catCi + CiCw) • (2-8) 

In this case, the one-particle wave function is the plane wave 



■ 2mTT(£ - 1) 
' N 



(2-9) 



where m is an integer that satisfies —N/2 + 1 < m < N/2. The corresponding 
one-particle energy is 

2m7r 

e^ = -2tcos^. (2-10) 
If is a multiple of four, the ground state energy at half filling is calculated as 

N/A 

r . = —2t cot 



E^= E ^m = -2tcot^. (2-11) 



m=-Af/4+l 

Thus, the finite-size correction to the energy per site 

Eff / 2 \ 2t TT 2t 2TTt , , 

-^-{ 1]= cot — + ^ 2-12) 

N \ TT J N N TT 3A^2 v; 

is of the order of 1/N'^. 

As verified in the above calculations, the finite-size correction to the energy per 
site E^ /N decreases faster for the system with the periodic boundary conditions 
than with the open boundary conditions. Regardless of this fact, the open boundary 
systems are often chosen in numerical studies by the density matrix renormalization 
group (DMRG) methocP^''^''^''^ because of the simplicity in numerical calcu- 
lation. It should be noted that for those systems that exhibits incommensurate 
modulation, the open boundary condition is more appropriate than the periodic 
boundary condition. Thus, it will be convenient if there is a way of decreasing the 
finite-size correction to E^ /N as fast as 1/A^^ also for the open boundary systems. 

§3. Spherical deformation 

We first consider the N-siie open boundary system described by the Hamiltonian 

N-l 

Hs = -tY^ sin (4c£+i + cl+i^q ■ (3-1) 



A^ 



Compared with the undeformed Hamiltonian Hq in Eq. (2-2), the strength of the 
hopping term is scaled by the factor = sin(^7r/A^), which decreases towards the 
system boundary as shown in Fig. 1. For a geometrical reason which we discuss in the 
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Fig. 1. A spherically deformed lattice, which contains {N = ll)-sites, drawn on the upper half 

of the circumference. Open circles denote lattice sites, where the angle of the £-th site is 
= k)''^/-^ ior i = 1, 2, . . . , N. The length of the vertical line shows the relative strength 
sin(£7r/iV) of the bond drawn by the thick arc between i-th and {£ + l)-th sites. 



next section, we call the modification from Hq to as the spherical deformation. 
We regard N, the number of sites on the upper half of the circumference shown in 
Fig. 1, as the system size. 

Let us observe the N dependence of the ground-state energy at half filling, 
where = {c\c^) = 1/2 is satisfied by the particle-hole symmetry. So far we have 
not obtained the analytic form of the one-particle wave function t/;^, except for the 
zero-energy state, and the corresponding one-particle eigenvalue for the deformed 
Hamiltonian H^. We therefore calculate them numerically by diagonalizing Hg in 
the one-particle subspace. We then obtain the expectation value {c\c^_^_l+c^^^-^^c^) and 
the ground state energy at half filling. In the following numerical calculations, 
we set t as the unit of the energy. 
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Fig. 2. The circles shows the expectation value (c]c£_|_i -|- c|_,_jq) of the spherically deformed lat- 
tice Fermion model defined by Hs when N = 1000. For comparison, we also plot the same 
expectation value for the undeformed case defined by Hq by the cross marks. 

Figure 2 shows (cjic^^^^ -|- c^^^-^c^) of the ground state when N = 1000. For com- 
parison, we also show the same quantity obtained by the undeformed Hamiltonian 
Hq of the same system size. As it is observed, the spherical deformation suppresses 
the position dependence in (c^^c^.^+c^^.^c^). In this sense we can say that the ground 
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state of is more uniform than that of Hq . 



One expects that the ground state energy Eg , which is the sum of negative 
one-particle eigenvalues 



N/2 



N-1 



N 

m=l e=i 
is nearly proportional to the sum of the bond strength 

N-l „ 

sin — = cot -—7 . 



(3-2) 



B 



N 



=1 



(3-3) 



It is also expected that the ratio /B^ rapidly converges to —2t/'K, which is the 
expectation value (cjc^,,.]^ + C£_,_j^c^) in the thermodynamic limit. Figure 3 shows 
E^/B^ and E^ /N with respect to 1/N'^. Obviously, the finit c-size corrections 
Ep /N + 2t/7r and E^ /B^ + 2t/TT are nearly proportional to l/N"^. In order to 
confirm this l/N"^ dependence, we show [E^ /N — Eg / B^]N'^ in Fig. 4, where the 
value converges to a constant in the limit N ^ 00. Calculating the ratio between 
\Ep /N + 2t/7r\ and |£'g'^ / B^ + 2t/7T\, we find that the former is twice as large as the 
latter in the limit N ^ 00. The result suggests that the spherically deformed A''-site 
system is related to a system of size 2N with periodic boundary conditions. 
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Fig. 3. The finite-size corrections to the energy per site at half filling n = 1/2. Crosses show 
/B'" and the open circles /N. 



We have considered the half-filled case. Away of the half filling we must include 
the chemical potential term, which is proportional to into the deformed Hamilto- 
nian. A natural way of introducing fj, is to put it into the bond operator h^^^^^, as 
stated in Eq. (1-2). From this extension we obtain the following Hamiltonian 




-tC^Cfij^l 



+ 4+i^e.+i 



(3-4) 
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Fig. 4. The convergence of [Ep /N - Es/B'^]N'^ with respect to 1/iV^. 
It is also possible to introduce the spherical deformation to the on-site terms as 

N-1 . N IP 

Hs = -t sm — [c'^c^+i + 4^iC^ J - 2^ sm — c'^c^ , (3-5) 

t=i e=i 

according to the height in Fig. 1 at each site. Note that both in Eq. (3-4) and 
in Eq. (3-5) give the same thermodynamic limit, and that the chemical potential 
terms do not commute with the kinetic energy in both cases. This is in contrast to 
the undeformed Hamiltonian 

N-1 N 

HQ = -tY^ (cjc^+i + 4+i^e) -f^Yl ' (3-6) 

e=i t=i 

where the chemical potential term is proportional to the total number of particles. 
Since we do not know the analytic formulation of one-particle energy of the deformed 
Hamiltonians in Eqs. (3-4) and (3-5), the relation between /x and the particle filling 
n = Ylei^^e)/-^ non-trivial. But we are interested in the cases where N is 
relatively large, therefore it is possible to use the relation n = —2t cos(7rn), which 
is satisfied by the undeformed Hamiltonian in Eq. (3-6) in the limit N 00, as a 
good approximation for fi for the spherically deformed system. 

Let us observe the occupation = {c^^Cf) with respect to the position i at 1/2, 
1/4 and 1/8 fillings, respectively, where the corresponding ^ is 0, — 2cos(7r/4), and 
— 2cos(7r/8). It is obvious that is always 1/2 at half filling, cquivalcntly, when 
/X = 0. Figure 5 shows calculated for in Eq. (3-4) when N = 1000. The 
particle distribution is almost uniform, since the ratio of the hopping strength and 
the chemical potential is independent on the position I on the lattice. Figure 6 shows 
near the boundary of the system. The oscillations in decay rapidly with the 
distance from the boundary. It should be noted that the amplitude of this small 
oscillation in the particle density decreases with increasing the system size N. 
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Fig. 



Position dependence in = {c\ci) near the system boundary at quarter filling. 



Figure 7 shows the finite-size correction to the energy per bond at 1/4 and 1/8 
filhngs, calculated for both in Eq. (3-4) and H'q^ in Eq. (3-5). As it is observed at 
half filling shown in Figs. 3 and 4, the correction is again proportional to l/N"^. We 
have thus confirmed the l/N"^ scaling for the correction to the ground-state energy 
per site of the spherically deformed lattice-free-Fermion model. 



§4. Geometrical interpretation 



There is a 2D classical system behind a ID quantum system, where the relation is 
called as the quantum-classical correspondence. We show that spherically deformed 
Hamiltonian i/g corresponds to a classical system on a sphere. We first consider the 
quantum-classical correspondence by way of the Trotter decomposition. For 
simplicity we consider the half- filled case (/i = 0) for the moment. Let us divide 
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Fig. 7. The finite-size corrections to the energy per site, where the crosses show /B^ calculated 
from the Hamiltonian Eq. (3.4), the open circles / from H'c. in Eq. (3.5). For com- 
parison we also show the correction E^ /N for the system with periodic boundary conditions 
by the triangles. 



in Eq. (3-4) into two parts 



(4-1) 



f=even 



=odd 



where we have used the notation h^^_^_i = —t (^cj'^^-fi + '^l+i'^y ^ where the de- 
formation factor is given by = sm^in/N). 

The imaginary time evolution of amount of /? is then expressed by the operator 
Q-l3Hg gy applying the Trotter decomposition to e~^^s ^ we obtain 



M 



(4-2) 



where M is the Trotter numbei'^''^ and Ap = /3/M. Looking at the structure of 
infinitesimal time evolution by 



exp ( -Z\/3 Ki+i J = exp ( - ^ (ApA^) h^^^^^ \ , (4-3) 

\ <?=even / \ £=even / 



we find that the quantity 

At^=A(3A^ (4-4) 

plays the role of the rescaled imaginary time. We can treat e~^^^2 in the same 
manner. It is possible to interpret Zir^ as a kind of proper tim^^ at the position 
i. Such interpretation leads us to an inhomogeneous time evolution on a (multiply 
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covered) sphere as shown in Fig. 8. This is the reason why we have used the term 
spherical deformation. Since the surface of the sphere is equivalent everywhere, it is 
natural to expect that the ground state of the spherically deformed Hamiltonian is 
approximately uniform. 




Fig. 8. Imaginary time evolution on a sphere. 



In the rest of this section we show the correspondence with the spherical geom- 
etry by taking the continuous limit to the lattice Hamiltonian or Hg. Consider 
a 1-particle state 



N 



(4-5) 



at time t. (Since we have been using the letter t for the hopping parameter, we use 
t for the time.) The real-time evolution of the wave function V'^(t) is described by 
the Schrodinger equation 



-t sin — V'^+i 



t sin 



{£ - l)7r 
N 



N 



(4-6) 



under the Hamiltonian H'^ in Eq. (3-5). Note that the difference between in 
Eq. (3-4) and in Eq. (3-5) is not relevant in the large N limit. There are two 
different continuous limits for this spatially discrete Schrodinger equation. We first 
consider the massive case where n is nearly equal to —2t. Introducing the notation 
= sin [{£ — i)7r/A^] , we can rewrite Eq. (4-6) by use of differentials 



Z^+i (V'^+i - i^e) - ft 



(4-7) 



where we have substituted the trivial relations V'^+i = (V'^+i ~ i^e) + V'^ i^e-i ~ 
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-{tp^ — + V'^. Using the relations 

1 



ft-l 



2 



2 



(4-8) 



we can further rewrite Eq. (4-7) as 



dt 



2 (■^^ 



2^^^ - l^fe'^e ■ 



(4-9) 



Now we introduce the lattice constant a = ttR/N, where R is the radius of 
the sphere. We also introduce the spacial co-ordinate x = a{i — which satisfies 
< X < ttR. Using these notations we rewrite V'^(t) as ^^{xjt), and as f{x) = 
sm(x/R) = sin 6, where 9 = x/R is the angle m,easured from the north pole. The 
continuous limit can be taken by increasing the number of sites keeping R constant, 
where the lattice constant a decreases with N. Simultaneously we increase the 
hopping parameter t so that the relation a^t = h? /2m always holds, where m is the 
particle mass, and h the Dirac constant. To prevent the divergence in the potential 
term, we adjust so that n+2t = —V is satisfied, where F is a finite constant. Using 
these parametrizations, we obtain the Schrodinger equation in continuous space 



/(^)^V'(aJ,t) 



+ Vf{x)i^{x,t). 



(4-10) 



This equation is derived from the Lagrangian 



d 

L = -ihip* [x, t) -Q^tp{x, t) + f{x) 



di;*ix,t)di;{x,t) 



2m dx 



dx 



+ ViIj*{x, t)tp{x,t) 



(4-11) 

where introduction of proper time r(a;,t) that satisfies dt = y^dr(a;,t) draws the 
following Lagrangian 



£ = /(x) 



in the x-t space. The action S is then written as 



,^ 9 , h'^ dip* dip _ ,' 
-ihip*—ip + — + ViP*iP 
OT 2m ox ox 



sin — drdx . 
R 



(4-12) 



(4-13) 



As it is seen, f{x) drdx = sin{x/R) drdx plays the role of the integral measure on the 
sphere of radius R. Note that the continuous limit for the field operator — >■ ip{x) 
can be taken in the same manner as in Eqs. (4-6)-(4-13) using the correspondence in 
Eq. (4-5). 
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Since f{x) = at the both ends, where x = and x = ttR, the continuous 
one-dimensional quantum system in Eqs. (4-10)-(4-13) does not effectively contain 
the system boundaries. We can observe the fact by way of the conformal mapping 

y = -R log cot (^) (4-14) 

from the sphere embedded in three dimensions onto the infinite plane. We have the 
relations 

sin — = 2 sin — cos — = ( cosh — ) (4-15) 

R 2R 2R \ RJ ^ ' 

and sm(x /R)dy = sinOdy = dx. The action <S on this infinite t-y plane is then 
written as 



sin ^ at / 2m \ sm ^ oy J \ sm oy 



sin^ 6* dtdy , 
(4-16) 

where sin 6 = (cosh ^) ^ is satisfied. The corresponding one-particle Hamiltonian 
is obtained as follows 

d ( . X d\ . X ^y d"^ 



2m dx 



sin— — -FF sin— = -— cosh ^—r-FF cosh . (4-17) 
\ RdxJ R 2m Rdy^ V RJ ^ ^ 

We can also formulate a massless limit, which appears in the case —2t < ji <2t 
where there is a Fermi surface, in the same manner as in Eqs. (4-6)-(4-9). In this case 
we substitute V'^ = e^**^^^^ to Eq. (4-6), where k and —k are, respectively, the Fermi 
wave number for the right and the left going modes. One finds that the quantity 
i' = at\s the leading order in the small lattice constant limit a — >^ 0, equivalently in 
the large N limit. Adjusting ji so that 2t cos k + ji = —V \s satisfied, we obtain the 
equation of motion 

ih—4>{x,t) = T'iivsuik—[f{x)(l){x,t)\ + Vf{x)(l){x,t), (4-18) 
for the continuous field 0(x,t). The corresponding Lagrangian in x-t plane is 

>C = f{x) 



d 0(j)* (x, t) 
ih(f)*{x,T)—(f){x,T) ± 211^ sink (f){x,T) ^— ^ h r)^(a;, r) 



(4-19) 

where we have used the fact that /(O) = f{2R) = 0. Similar to Equations (4-13)- 
(4-16), we can consider the conformal mapping for this massless case. The l/N"^ 
dependence of the corrections to the ground-state energy per site might be explained 
by the boundary conformal field theory, where we leave the conjectures for the future 
study. 

§5. Conclusions and discussions 



We have investigated the ground state of the spherically deformed ID free 
Fermion system, for both at the half filling and away of the half filling. The finite- 
size correction to the energy per site is of the order of l/AT^ for both cases. The 
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reason for such fast convergence is qualitatively explained by the quantum-classical 
correspondence, where the spherically deformed Hamiltonians essentially correspond 
to classical fields on a sphere. In such a sense the spherically deformed system does 
not contain the system boundary. 

Interest in the spherical deformation rests in dynamical properties. We con- 
jecture that a moving one-particle wave packet on the spherically deformed lattice 
oscillates nearly harmonically as a consequence of the circulation on the sphere. 
The oscillation may be also explained by a continuous refraction caused by a slower 
dynamics near the both ends of the system. 

As a generalizations of the spherically deformed Hamiltonian Hg in Eq. (1-7), 
one can consider a decoupled Hamiltonian 

Hsin = sm— I + ) , (5-1) 

e=i ^ ^ 

where £ = 2N -|- 1 is equivalent to ^ = 1, for a system of size 2N. The differential 
with respect to £ draws 

Hcos = cos— I + ^ I , (5-2) 

which is again the decoupled Hamiltonian when is an even number.!^ Both ^^jj^ 
and -ffcos seems to be generators of rotation on a kind of discrete sphere. Their 
commutation relation would be discussed elsewhere. 

If one is interested in the estimation of the excitation gap, the spherical defor- 
mation is not appropriate. This is because weak bonds near the system boundary 
induce spurious low-energy excitations. For this purpose, the hyperbolic deformation 
is more appropriate.'2D''E3 The quantum-classical correspondence discussed in this 
article can be also considered for the hyperbolic deformation, which would deduce 
continuous field model on the Poincare disc. 
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